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Abstract 

The inclusive production of £)* ± (2010) mesons in deep-inelastic scattering is studied with 
the HI detector at HERA. In the kinematic region 1 < Q 2 < 100 GeV 2 and 0.05 < y < 0.7 
an e + p cross section for inclusive D meson production of 8.50±0.42 (statOj^oo (syst) nb 
is measured in the visible range p t D* > 1-5 GeV and |r/_o* \ < 1-5. Single and double 
differential inclusive D* ± meson cross sections are compared to perturbative QCD calcu- 
lations in two different evolution schemes. The charm contribution to the proton structure, 
F% (x,Q 2 ), is determined by extrapolating the visible charm cross section to the full phase 
space. This contribution is found to rise from about 10% at Q 2 = 1.5 GeV 2 to more than 
25% at Q 2 = 60 GeV 2 corresponding to x values ranging from 5 • 10~ 5 to 3 • 10" 3 . 
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1 Introduction 



Results on inclusive D*^ meson production in deep-inelastic ep scattering (DIS) and on the 
charm contribution to the proton structure function, F|, at HERA have been published by the 
HI and the ZEUS collaborations [JT] — 13|] - These data, together with earlier fixed target data [Q], 
have shown clear evidence that the dynamics of charm production in ep scattering is described 
by the photon gluon fusion process, which is sensitive to the gluon density in the proton [|5|] and 
allows its universality to be tested. 

Early results on F 2 C from the HI experiment [[]]] were based on an integrated luminosity of 
3 pb 1 collected during the 1994 HERA running and were therefore statistically limited. The 
current analysis uses data from the 1996 and 1997 HERA running periods, yielding a signifi- 
cantly larger integrated luminosity of 18.6 pb -1 . Furthermore, the improved instrumentation in 
the backward region of the HI detector enables the kinematic range in four-momentum transfer 
squared to the virtual photon, Q 2 , to be significantly extended down to 1 GeV 2 . Hence, more 
precise tests of perturbative QCD (pQCD) become possible. 

This paper is organized as follows: a short discussion of the different approaches to open 
charm production in perturbative QCD calculations is followed by a description of the experi- 
mental set-up and details of the analysis; the inclusive cross sections for D* ± meson production 
are then presented and compared to QCD predictions. Finally, they are used to derive the charm 
contribution to the proton structure function, F 2 C . 



2 Models of Open Charm Production 

The description of open heavy flavour production in electron proton collisions is based on per- 
turbative QCD. In leading order (LO), the photon gluon fusion process (jg — > QQ) is the 
dominant contribution [[[]]. Next-to-leading order (NLO) calculations in several schemes are 
available [|6|-|T0|]. All approaches assume that Q 2 and the heavy quark mass itlq provide a hard 
enough scale to allow the applicability of pQCD and to guarantee the validity of the factoriza- 
tion theorem. 

Here, the "massive approach" is adopted, i.e. a fixed order calculation with massive quarks 
assuming three active flavours in the proton. The momentum densities of the three light quarks 
and the gluon in the proton are evolved by the DGLAP equation [JTT|] . The heavy quarks are 
assumed to be produced only at the perturbative level ^\ via photon gluon fusion. These cal- 
culations are considered reliable in the regime Q 2 m tuq. However, they break down at some 
scale Q 2 3> tuq due to large logarithms ~ ln(Q 2 / 'm 2 ^) . 

Based on fixed order a 2 calculations in the coefficient functions [j^] programs for different 
applications were developed. The Riemersma et al. program [0] can be used to calculate inclu- 
sive quantities of heavy quark production, like F£(x, Q 2 ), while the HVQDIS program |||[T2|] 



allows the calculation of exclusive quantities by providing the four-momenta of the outgoing 
partons. In the version of the program used here charmed quarks are fragmented in the pho- 
ton - proton centre of mass frame into Z)* ± mesons using the Peterson fragmentation func- 
tion [p"3"|], which is controlled by a single parameter e c . In addition, to account for the experi- 
mentally observed p t smearing of hadrons with respect to the quark direction, the meson 
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has been given a transverse momentum p t with respect to the charm quark, according to the 
function p t ■ exp(—ap t ). The parameter a is chosen such that an average transverse momentum 
(p t ) ~ 350 MeV is obtained as observed in e + e~ data [|T3]]. With this procedure it becomes pos- 



sible to calculate differential inclusive D* ± meson cross sections in the experimentally visible 
phase space region. 

The CCFM evolution equation [ |T5| ] is expected to be more appropriate to describe the parton 
evolution at small x. In the parton cascade, gluons are emitted in an angular ordered manner 
to account for coherence effects. Due to this angular ordering, the gluon distribution depends 
on the maximum allowed angle in addition to the momentum fraction x and the transverse 
momentum of the propagator gluon. The cross section is then calculated according to the k t - 
factorization theorem by convoluting the unintegrated gluon density with the off-shell photon 
gluon fusion matrix element with massive quarks for the hard scattering process. 

It has been shown previously [ |16| ] that F 2 and F 2 C can be reasonably well described within the 
CCFM framework. In addition a solution of the CCFM equation has been obtained recently [ [T7| ] 
from a fit to F 2 which is able to describe the cross section for forward jet production, where sig- 
nificant differences to the expectation in the DGLAP evolution scheme are seen. Using this 
solution the hadron level Monte Carlo generator CASCADE has been developed [|18p. This 
allows the full generation of charm events including the initial state gluon radiation according 
to the CCFM equation and the fragmentation of partons by the Lund String model. The frag- 
mentation of charmed quarks to D* ± mesons is performed using the Peterson fragmentation 
function. 



3 Detector and Simulation 



The data have been collected with the HI detector [ |T9| ] at HERA during the running periods of 
1996 and 1997 when HERA operated with 27.5 GeV positrons and 820 GeV protons colliding 
at a centre of mass energy of y/s = 300 GeV. The following detector components are important 
for this analysis. The scattered positron is identified and measured in the SpaCal POTJ, a lead- 
scintillating fibre calorimeter situated in the backward regionQof the HI detector. The SpaCal 
also provides time- of- flight information for trigger purposes. A four double-layer backward 



drift chamber (BDC) Q21| ] is mounted in front of the SpaCal in order to improve the angular 



measurement of the scattered positron. Charged particle tracks are reconstructed by two cylin- 



drical central jet drift chambers (CJC) [|TJ, f22p placed concentrically around the beam-line in 
a homogeneous magnetic field of 1.15 Tesla. The CJC also provides trigger information [ ]23| ] 
based on the detection of track segments. Double layers of cylindrical multi-wire proportional 
chambers (MWPC) [ J24] ] for triggering purposes are positioned inside and in-between the two jet 
chambers. The luminosity is determined from the rate of the Bethe-Heitler reaction ep — > ep'j. 

Monte Carlo simulation programs are used to simulate detector effects and to estimate the 
systematic uncertainties associated with the measurement. For the determination of the accep- 
tance of the detector and the D*^ selection efficiencies, heavy flavour (charm and bottom) DIS 



1 The positive z-axis of the HI reference frame, which defines the forward direction, is given by the outgoing 
proton direction. 
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events are generated using the AROMA 2.2 p5| ] program. This program, which is based on 
the DGLAP evolution scheme, simulates neutral current heavy quark production via photon 
gluon fusion in leading order QCD including parton showers and heavy quark mass effects. 
The mass of the charm quark is chosen to be m c = 1.5 GeV while the factorization and renor- 
malization scales are set to \i = V~§, where s denotes the square of the invariant mass of the 
heavy quark system. The GRV94-LO [ p6| ] parton density functions (PDF's) are used for the 
proton. Hadronization is performed in the Lund String Model [fZ7p, as implemented in JETSET 
7.4 [|28|]. The momentum fraction of the charm quark carried by the D* ± meson is determined 
according to the Peterson model [JT3p with the fragmentation parameter e c = 0.078 []Z9[]. The 
influence of the details of the fragmentation process on the acceptances and efficiencies has 
been investigated by (a) varying the Peterson fragmentation parameter between e c = 0.035, as 
favoured in [|30|], and e c = 0.1 which seems to yield a better description of the hadronic final 
state in D*^ 1 events, (b) applying the symmetric Lund fragmentation function $T/\\ also to the 
D* ± mesons and (c) using the HERWIG [ |JT] ] program which is based on the cluster hadroniza- 
tion model [J32[]. The inaccuracy due to the uncertainty in the QCD parameters is studied by 
varying the charm quark mass m c and by changing the factorization and renormalization scales 
to /i = \/Q 2 + 4m;?. The dependence of the acceptances and efficiencies on the QCD evolu- 
tion scheme has been determined also by using the CASCADE [ |17| J event generator. Finally, 
the influence of QED radiation on the efficiency is determined using the RAPGAP [[33]] pro- 
gram interfaced to HERACLES 4.1 [|34j]. All Monte Carlo generated events are fed into the 
GEANT p5| ] based simulation of the HI detector and are subjected to the same reconstruction 
and analysis chain as used for the data. 



4 Kinematics 

This analysis is restricted to those DIS events which have a scattered positron detected in the 
backward region of the detector. At fixed center of mass energy y/s the kinematics of the inclu- 
sive scattering process ep — ► eX can be completely determined by any two of the independent 
Lorentz invariant variables: the Bjorken scaling variable x, the lepton inelasticity y, the four- 
momentum squared Q 2 = —q 2 of the virtual photon and the invariant mass squared W 2 of the 
hadronic final state. In this analysis, these variables are determined from the measurement of 
the energy E' e and the polar angle e of the scattered positron according to the expressions 

Q 2 = AE e E' e cos 2 (^) 

Q 2 
x = — 

ys 

where s = AE e E p and E e and E p denote the energies of the incoming positron and proton, 
respectively (the positron and proton masses are neglected). 



y = 1 sm — 



(1) 



W 2 = Q' 



1 — x 



x 
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5 Event Selection 



The events for this analysis were triggered by a coincidence of an electromagnetic cluster in the 
SpaCal with a charged track signal from the CJC and a vertex which is coarsely reconstructed 
from the MWPC information. The positron is identified as the most energetic cluster with 
E' e > 8 GeV as described in [|45|]. The cluster radius^ is required to be less than 4 cm, consistent 
with an electromagnetic energy deposition, and the cluster center of gravity is required to be 
within 1 .5 cm of the extrapolation of a charged track segment from the backward drift chamber 
BDC. The geometrical acceptance of the SpaCal and BDC imposes a limitation on the positron 
scattering angle of e < 177.5°. These limits and requirements restrict the accessible range in 
Q 2 from 1 GeV 2 to 100 GeV 2 and in the lepton inelasticity to y < 0.7. To measure the event 
kinematic quantities with sufficiently good resolution y is further constrained to y > 0.05. Good 
agreement is observed for all quantities related to the scattered positron between data and the 
prediction of the AROMA Monte Carlo simulation. The contribution due to photoproduction 
background, i.e. Q 2 < 1 GeV 2 , is everywhere smaller than 1% in the selected kinematic region. 

Charm production is identified by the reconstruction of mesons in the decay chain 

D* + - D\+ ow - {K-k+)k+ ow (+c.c.) (2) 

using the D* — D° mass difference method [TJ^]. The decay products are detected in the central 
track detector. For each accepted track particle identification is applied using the measurement 
of the energy loss, dE/dx, in the central track detector. In order to reconstruct a D*^ 1 meson 
candidate, unlike-sign charged tracks are first combined to form i^ T 7r ± pairs in which one of the 
particles should be consistent with a kaon and the other with a pion according to their dE/dx 
measurements. Among all possible oppositely charged i^ =F 7r ± pairs, those with an invariant 
mass consistent within ±70 MeV of the D° mass are combined with a track of a second pion 
candidate ("vr^ OJj; ") having a charge opposite in sign to that of the kaon. In Fig. [T| a clear peak 
is observed in the distribution of the mass difference Am = m^?™ — m Kn around the nominal 
£)*± _ £)0 mass difference of 145.4 MeV. A fit to this distribution using a Gaussian for the 
signal and a term (Am — m n ) a for the background yields a total of 973 ± 40 D*^ 1 mesons in the 
acceptance range of pseudorapidity^ \rj K7T7T \ < 1.5 and transverse momentum pud™ > 1-5 GeV. 



6 Inclusive Cross Sections 

The integrated and differential Born level cross sections for meson production in DIS are 
calculated from the observed number N D *± of candidates, according to 

-.e+ir**)- N °- ±( \Z\ (3) 

Lint ■ D ■ € ■ [l + rad J 

Here, r stands for the contribution of reflections in the D° mass window, coming from D° 
channels other than the one studied in this analysis. The value of r amounts on average to 

2 The cluster radius is defined as log(Ei) ■ di/ J2i 1°SC^») where the sum runs over all cells in the cluster: 
Ei is the normalized energy of the cell i and di is the distance of the cell i from the cluster centre of gravity. 
3 The pseudorapidity of a particle is defined as rj = — lntan(9/2). 
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about 0.03. The integrated luminosity is denoted by C int while B refers to the branching ratio 
B = B(D* + -> D°n+)-B(D° -> A^vr + ) = 0.0259±0.0006 ©. The detection efficiency e is 
estimated to be 22.5% using AROMA. The radiative correction <5 md which correct to the single 
photon exchange cross sections are obtained from the program HECTOR [|3^]. Depending on 
the kinematic region 5 rad varies from +0.11 at small x and Q 2 to -0.02 at large Q 2 . For the 
integrated visible cross section it averages to 0.03. 



6.1 Integrated Cross Section 

The inclusive cross section for meson production in the kinematic region 1 < Q 2 < 
100 GeV 2 and 0.05 < y < 0.7, and in the visible D*^ range \ij D *\ < 1.5 and p tD * > 1.5 GeV 
is found to be 

avis(e + p -> e + D* ± X) = 8.50 ± 0.42(stat.)t^ (syst.) ± 0.65(model) nb. 

The errors refer to those from statistics, experimental systematics and additional systematics 
related to the changes in efficiency obtained by using different Monte Carlo generators and 
varying the model parameters. 

The experimental systematic uncertainties are summarized in Table |T]. The largest contri- 
bution is due to the uncertainty in the track reconstruction efficiency. Other important sources 
include uncertainties in the extraction of the D*^ signal, i.e. the determination of the back- 
ground shape in the Am distribution and the D° mass resolution. The uncertainties due to 
model dependencies, as summarized in Table |2|, include the incomplete understanding of the 
fragmentation process, the uncertainty due to the charm quark mass, the sensitivity to the fac- 
torization and renormalization scales and the change of acceptance due to QED effects at the 
positron vertex. The largest effect on the efficiency is observed by changing the charm quark 
mass from m c = 1.5 GeV in the reference Monte Carlo dataset to m c = 1.3 GeV and by 
changing the fragmentation models and their parameters. 

The visible inclusive meson production cross section has been calculated in the NLO 
DGLAP scheme with the HVQDIS program using the GRV98-HO parton densities in the pro- 



ton [|39|]. The predictions range from 5.17 nb for a charm quark mass m c = 1.5 GeV and 



Peterson fragmentation parameter e c = 0.10 to 7.02 nb for m c = 1.3 GeV and e c = 0.035. 



The hadronization fraction f(c D* + ) = 0.233 ± 0.010 ± 0.011 [gOJ has been used. For 
the same variation of m c and e c , calculations based on the CCFM evolution, as implemented 
in the CASCADE program, yield a significantly higher cross section of 8.04 nb and 10.77 nb, 
respectively. 

Disregarding the small differences in the kinematic range, good agreement is observed in the 
inclusive meson production cross section with the result obtained by the ZEUS experiment 
[0]. The measured value of this cross section agrees better with the CASCADE prediction than 
with that from HVQDIS. In previous publications [|l],|5|] HI reported much better agreement 
between data and predictions from the HVQDIS program. The larger difference obtained now 
is due to the new determination of the charm quark hadronization fraction f(c — > D* + ) which 
is 16% smaller than the previous value. 
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6.2 Differential Cross Sections 



In Fig. [Z]the inclusive single differential D*^ 1 cross sections in the visible region are shown as 
a function of the event variables W, x and Q 2 and as a function of the D*^ observables ptD*, 
j] D * and the inelasticity z D * = P ■ p D *j 1 P ■ q = (E — p z ) D * /2yE e , where P, q and p D * denote 
the four-momenta of the incoming proton, the exchanged photon and the observed _D* ± meson, 
respectively. A bin by bin correction to account for QED radiation has been applied. 

Fig. ^ also includes the expectations from the HVQDIS program using the GRV98-HO 
parton density parameterization. The renormalization scale and the factorization scale are set 
to [i = a/<5 2 + 4m 2 .. The charm quark mass and the fragmentation parameter have been varied 
from m c = 1.3 GeV and e c = 0.035 to m c = 1.5 GeV and e c = 0.10. The dark shaded band 
indicates the uncertainties in the predictions due to these variations. Although the predicted 
visible cross section is smaller than experimentally observed, the agreement with the data in the 
shapes of the different single differential cross sections is reasonable. A significant difference 
is observed in the da/dr/ cross section. For > the measured D*^ meson production cross 
section is larger than predicted by the calculation. Since in the boson gluon fusion process the 
forward region (i] D * > 0) is correlated with small Zr>* a similar discrepancy between data and 
theory is observed at small Zd*- 

A possible cause of this deviation could be the simplified grafting of fragmentation onto 
the HVQDIS program. This approach does not account for the colour force between the charm 
quark and the proton remnant which is expected to result in a drag of the D*^ meson from the 
original charm quark direction towards the proton direction. To quantify this 'beam drag effect' 
PH], a mapping function from the r] c -p tc space to the T]D*-ptD* space has been constructed 
using the AROMA Monte Carlo program which includes such effects. This function has then 
been used instead of the Peterson fragmentation with transverse momentum smearing in the 
HVQDIS program. No significant change in the rjo* and zd* distributions has been observed 
by this procedure compared to our original treatment of fragmentation. A better description 
of the r/D* distribution is obtained, however, when using the HERWIG program to extract the 
mapping function, at the expense of a 10-15% reduction in the visible cross section prediction. 
It is therefore concluded that the differences between the measurements and the predictions 
from the HVQDIS program can not be explained by the absence of colour drag effects in these 
calculations. 

Fig. [2] also presents the predictions of the CASCADE program with the same variations of 
the charm quark mass and the fragmentation parameter. The expectations from the CASCADE 
program are found to agree better with the data in general and especially in the positive r] region. 

In order to enable the study of correlations among the observables in D* ± meson production, 
Figs. |3| and |] show the double differential inclusive cross sections. It is evident that the 
excess observed in the data with respect to the HVQDIS expectation at large pseudorapidities 
(0.5 < t]d* < 1-5) is independent of Q 2 and is concentrated at small p t r>* and small zd*. It is 
especially in this phase space region where the CASCADE program better represents the data. 
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7 Charm Contribution to the Proton Structure Function 



The charm contribution, F 2 c (x, Q 2 ), to the proton structure function is obtained by using the 
expression for the one photon exchange cross section for charm production 

^=w (1+(1 - y)2)F ^ Q2) ' (4) 

where the contribution of the longitudinal structure function is neglected. The visible inclusive 
D* ± cross sections a^(x,Q 2 ) in bins of x and Q 2 are converted to a bin center corrected 
F 2 6XP ((^), (Q 2 )) by the relation: 



exp/ /~)2\ 

^ cexp «*>, (Q 2 )) = a Zo * J, ■ FS theo ((x), (Q 2 )) , (5) 

^vis ( x ,Q ) 



where and F| theo are the theoretical predictions from the model under consideration. The 
measured values of the visible cross sections, o\)* p (x, Q 2 ), are listed in Table |3|. Following the 
same line as in previous publications the HVQDIS program by Harris and Smith JT2| ] 

and the program of Riemersma et al. [^p are used to calculate these quantities in the NLO 
DGLAP scheme. In the kinematic range of the current analysis the beauty contribution to the 
proton structure, F 2 b , is expected to be of the order of 1 to 2% of F 2 C [ft2"|,ft3|] and is therefore 
neglected Q. 

In Fig. |]a F 2 C is shown as a function of x for different values of Q 2 as extracted from the 
inclusive D*^ cross sections using m c = 1.4 GeV. The systematic error on the data points 
includes those described in Sec. |(T2| as well as additional errors coming from the extrapolation 
in Equation |[ The bands show the predictions based on the gluon density extracted by the HI 
NLO DGLAP fit to the inclusive F 2 measurement 0431]. The width of each band reflects the 



total uncertainty of the prediction resulting from the uncertainties on this fit, thereby exploiting 
the full correlations arising from the constraints of the inclusive F 2 measurement. The influ- 



ence of all the individual sources of uncertainties considered in Ref. [|45j] for the determination 
of the gluon density and the strong coupling constant a s has also been investigated here. The 
most relevant variations with respect to theoretical calculations are the variation of the strong 
coupling constant a s in the range 0.113 < a s < 0.167, of the factorization and renormalization 
scale [i in the range 0.5 • (Q 2 + Am 2 ) < // < 2 • (Q 2 + Ami) an( ^ °f tne charm quark mass in 
the range 1.3 < m c < 1.5 GeV. The dominant sources of uncertainties are the experimental 
error on the F 2 measurement at very small x and the insufficient knowledge of the charm quark 
mass in the range of the direct F| measurements. For the displayed bands the different contri- 
butions are added in quadrature. Fig. |5Ja also includes the results of [TJ] for comparable values 
of Q 2 . These measurements suggest a steeper rise at small Q 2 towards small x than expected 
from the calculations based on the gluon density in the proton extracted from the inclusive F 2 
measurement. 

The extraction of F 2 according to Equation |5J is faced with an intrinsic problem. The mea- 
surement covers about 30% of the total phase space for charm production and the estimation 

4 If, however, the beauty cross section turns out to be large [p4[], its contribution to _D* ± meson production may 
have to be subtracted from the visible inclusive cross section prior to the determination of F£. 
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of this acceptance fraction depends significantly on the underlying model. To be more explicit, 
two different calculations may yield the same value for F| theo ( (x) , (Q 2 ) ) but may have different 
acceptances. They may then predict for a bin in x and Q 2 different cross sections a^°(x, Q 2 ) 
with consequent different values for F 2 cxp . To investigate the model dependence F 2 cxp is also 
determined using the CASCADE program with m c = 1.4 GeV and the results are shown in 
Fig. |]b. The figure includes also the prediction according to the CCFM evolution. Here the 
bands indicate the uncertainty on this prediction due to the variation of the charm quark mass. 
The comparison of Figs. |]a and reveals a steeper rise in the predicted charm contribution 
to the proton structure function at small x in the CCFM evolution than obtained by the NLO 
DGLAP evolution. Using the acceptances and efficiencies calculated from the CASCADE pro- 
gram the measured values of F 2 are found to be systematically smaller than those determined 
with the HVQDIS program. The largest differences (up to ~ 20%) are observed at small x 
values. 

In Fig. |6| F 2 is shown as a function of Q 2 for different values of x using the acceptances 
as calculated with HVQDIS. As in Fig. |]the bands indicate the full uncertainty in the DGLAP 
NLO predictions for a central value of the charm quark mass of 1.4 GeV using the gluon density 
extracted from the fit to the inclusive F 2 measurement. The full line shows the DGLAP NLO 
prediction using the gluon density from GRV98-HO. Taking into account the different data sets 
used for the determinations of the gluon densities, the agreement of the different calculations 
is reasonable. The data show a steep rise of F 2 with Q 2 . The slope, <9F 2 /<91n Q 2 , contributes 
roughly half of the slope of the inclusive structure function, dF 2 /dln Q 2 , measured at the cen- 
tral Q 2 of each x bin in [ p3Q in the range 0.0002 < x < 0.002. This steep rise is reasonably 
well reproduced by the NLO DGLAP calculations. 

In Fig. [7] the ratio of F 2 to the inclusive F 2 [@3|] is shown as a function of x for different 
values of Q 2 . The contribution of charm production to the total F 2 rises from about 10 % at 
Q 2 = 1.5 GeV 2 and x « 10~ 4 to more than 25 % at Q 2 > 25 GeV 2 and x > 5 • 10~ 4 . Although 
this behaviour agrees with expectation, at small x the measured ratio of F 2 /F 2 is larger than 
predicted in the NLO DGLAP scheme. 



8 Conclusions 

New measurements of differential cross sections for inclusive -D* ± production in deep-inelastic 
ep scattering are presented. These are compared with predictions based on both NLO DGLAP 
and CCFM formalisms, the former made using the HVQDIS program [ ]T2"| ] and the latter using 
the CASCADE model [|T7p, The predictions made using DGLAP formalism tend to undershoot 
the data, particularly for small D* transverse momenta, ptD*, and positive D* pseudorapidities, 
i] D * . The expectations of the CCFM based model are in better agreement with the data. 

Extrapolation of the visible D* cross section to the full p t n* and t]d* phase space allows 
extraction of F 2 , the contribution of charm to the proton structure function F 2 . These extrapo- 
lations are seen to depend on the formalism used: that based on the NLO DGLAP formalism 
typically produces a larger result for F| than that made using the CCFM approach. Both re- 
sults are presented to allow consistent comparisons using either formalism. The kinematic 
range presented has been extended to lower Q 2 than shown in the previous HI study, namely 
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Q 2 = 1 GeV 2 , and thereby to lower x. The F 2 C measurements show large scaling violations and 
a steep rise of F% with decreasing x. This rise tends to be steeper than expected from the NLO 
DGLAP calculations, but agrees well with the CCFM based expectations. Both approaches 
for the extraction of F| show that the contribution of charm production to F 2 exceeds 25% for 
Q 2 >25GeV 2 . 
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Experimental Systematic Uncertainties 



Trigger efficiency 




±0.02 


Track detector efficiency 


+0.075 


-0.034 


dE/dx measurement 




±0.03 


Pt(7T s low)-CUt 




+0.025 


Background shape 




±0.05 


D° Mass resolution 


+0.04 


-0.025 


Reflections 




±0.015 


Event kinematics 




±0.04 


Luminosity measurement 




±0.015 


7P contribution 




<0.004 


Branching ratio 




± 0.025 




+0.12 


-0.09 



Table 1: Summary of the fractional experimental systematic uncertainties of the inclusive D* 
meson cross section. 



Model Uncertainties 

Fragmentation model +0.035 —0.07 

Charm quark mass +0.07 —0.02 

Scale fi = y/Q 2 + 4m 2 c -0.015 

QED radiation -0.025 

±0.08 



Table 2: Summary of the fractional model dependent uncertainties of the inclusive meson 
cross sections. 
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l0g G?V 2 


(Q 2 ) 


log (; 




(x) 


CTvis 




Ft 




Ft 






[GeV 2 ] 






[io- 3 ] 


[nb] 




(DGLAP) 


(CCFM) 





1.5 


-4.8 to 


-4.2 


0.05 


1.11 ±0.19 


+0.15 
-0.28 


0.114 ±0.019 


+0.019 
-0.030 


0.094 ±0.016 


+0.014 
-0.024 


to 0.375 




-4.2 to 


-3.0 


0.20 


1.44 ±0.28 


+0.55 
-0.34 


0.053 ±0.010 


+0.021 
-0.013 


0.056 ±0.011 


+0.022 
-0.013 


0.375 


3.5 


-4.6 to 


-3.6 


0.13 


0.95 ±0.13 


+0.14 
-0.13 


0.161 ±0.022 


+0.023 
-0.027 


0.159 ±0.021 


+0.024 
-0.022 


to 0.625 




-3.6 to 


-3.0 


0.50 


0.41 ±0.07 


+0.11 
-0.09 


0.082 ±0.015 


+0.022 
-0.018 


0.082 ±0.015 


+0.022 
-0.018 


0.625 


6.5 


-4.0 to 


-3.6 


0.20 


0.72 ±0.10 


+0.27 
-0.08 


0.309 ±0.041 


+0.116 
-0.042 


0.256 ±0.034 


+0.097 
-0.031 


to 1 




-3.6 to 


-3.0 


0.50 


1.02 ±0.11 


+0.21 
-0.10 


0.151 ±0.017 


+0.032 
-0.015 


0.153 ±0.017 


+0.032 
-0.015 


1 


12 


-3.6 to 


-3.0 


0.50 


0.59 ±0.09 


+0.11 
-0.11 


0.276 ±0.040 


+0.052 
-0.052 


0.253 ±0.037 


+0.048 
-0.047 


to 1.25 




-3.0 to 


-2.0 


2.00 


0.43 ±0.07 


+0.16 
-0.04 


0.174 ±0.029 


+0.064 
-0.018 


0.194 ±0.032 


+0.071 
-0.021 


1.25 


25 


-3.6 to 


-3.0 


0.50 


0.33 ±0.06 


+0.08 
-0.06 


0.508 ±0.096 


+0.125 
-0.087 


0.422 ±0.079 


+0.105 
-0.071 


to 1.5 




-3.0 to 


-2.0 


2.00 


0.54 ±0.07 


+0.07 
-0.09 


0.278 ±0.037 


+0.035 
-0.048 


0.273 ±0.036 


+0.036 
-0.047 


1.5 


60 


-3.0 to 


-2.6 


2.00 


0.31 ±0.06 


+0.07 
-0.03 


0.348 ±0.066 


+0.077 
-0.034 


0.362 ±0.068 


+0.081 
-0.035 


to 2 




-2.6 to 


-2.0 


3.16 


0.35 ±0.06 


+0.09 
-0.06 


0.272 ±0.050 


+0.067 
-0.044 


0.270 ±0.050 


+0.069 
-0.043 



Table 3: Inclusive D* ± cross section a vis in bins in x and Q 2 for the visible range \rj D * | < 1.5 
and p t D* > 1-5 GeV as extracted in the DGLAP scheme. The values given in the columns 
denoted by (Q 2 ) and (x) are the bin centres at which the values of F| are given. The values of 
F| extracted in both the DGLAP and CCFM schemes are also presented. 



16 



300 



> 

M 

o 
o 

IT) 

QJ 
Oh 

8 200 

T3 



100 








HI 



□ ( J ftT ± 7T ± )7rf 




0.14 



0.15 



0.16 



m Kmz~ m Kn 



0.17 



[GeV] 
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Figure 3: Double differential inclusive cross section d 2 a/dr]dQ 2 and d 2 a/dp t dQ 2 in bins 
r}£>* and p t d* (see figure for details). 
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Figure 5: F 2 C , as derived from the inclusive D*^ meson analysis (a) in the framework of 
NLO DGLAP and (b) in the framework of CCFM, both for m c = 1.4 GeV. The error bars 
on the HI data points refer to the statistical (inner) and the total (outer) error, respectively. 
In (a) the shaded bands represent the predictions of F 2 C from the HI NLO DGLAP fit to the 
inclusive F 2 measurements including all the uncertainties described in the text. The dominant 
contribution arises from the uncertainty of m c . The ZEUS measurements [|3p are shown for 
comparable values of Q 2 indicated in parantheses (see [^] for a discussion of the extrapolation 
uncertainties). In (b) the bands represent the expectation of F 2 C from the fit to the inclusive F 2 
in the CCFM scheme including only the variation of m c which in both cases ranges from 1.3 to 
1.5 GeV. 
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Figure 6: F 2 C as derived from the inclusive D*^ 1 meson production cross section as a function 
of Q 2 for different values of x. The error bars refer to the statistical (inner) and the total (outer) 
error, respectively. The shaded bands represent the predictions of the NLO DGLAP evolution 
based on the parton densities in the proton obtained by the fit to the inclusive F 2 for m c = 
1.4 GeV including all uncertainties described in the text. The black lines show the predictions 
in the NLO DGLAP scheme for m c = 1.4 GeV using the gluon density from GRV98-HO. 
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Figure 7: The ratio of F 2 C over F 2 as derived from the inclusive D*^ meson analysis as a function 
of x for different values of Q 2 . The error bars refer to the statistical (inner) and the total (outer) 
error, respectively. The shaded bands represent the predictions of the NLO DGLAP evolution 
based on the parton densities in the proton obtained by the fit to the inclusive F 2 for a central 
charm quark mass of 1.4 GeV including all uncertainties. 
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